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Abstract 

In the BMN approach to A/" = 4 SYM a large class of correlators of interest are express- 
ible in terms of expectation values of traces of words in a zero- dimensional Gaussian 
complex matrix model. We develop a loop-equation based, analytic strategy for eval- 
uating such expectation values to any order in the genus expansion. We reproduce 
the expectation values which were needed for the calculation of the one-loop, genus 
one correction to the anomalous dimension of BMN-operators and which were earlier 
obtained by combinatorial means. Furthermore, we present the expectation values 
needed for the calculation of the one-loop, genus two correction. 
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1 Introduction 



Recent progress in string and gauge theory has brought to hght an interesting 

pp-wave/BMN-correspondence which is a special version of the celebrated AdS/CFT 
correspondence. The pp-wave is a ten-dimensional geometry which can be obtained as 
a Penrose limit of AdSs x S5 and which constitutes a background where it is possible to 
quantize type IIB string theory in light cone gauge . BMN stands for Berenstein, 

Maldacena and Nastase who identified the gauge theory dual as a special sector of 
M = 4 SYM based on gauge group SU(N) with a certain limit understood. 

The BMN sector of A/" = 4 SYM consists of operators which carry a large R-charge, 
J, associated with a selected SO (2) sub-group of the full SO (6) R-symmetry group and 
for which A — J is finite where A is the conformal dimension. It has been argued that 
the quantum corrections to correlation functions involving such operators only depend 
on qym via the parameter A' = {gyM^)/ J"^ ^"^^ the BMN limit is given by the 

scaling prescription 

Qym fixed, J, — > 00 with (72 = — fixed (1.1) 

which in particular renders A' finite.0 As shown in [^H, despite being a large-N 
limit the BMN limit is not a planar limit. Diagrams of all genera survive the limit 
and contributions of genus h are weighted by a factor {g2Y^- One could say that the 
BMN approach to A/" = 4 SYM introduces a new 't Hooft expansion with a new gauge 
coupling constant A' and a new genus counting parameter, g2. However, the BMN limit 
is not a new 't Hooft limit because the genus counting parameter remains finite as the 
limit is taken. Rather, the BMN limit is an interesting new double scaling limit much 



like the one encountered in the study of 2D quantum gravity |10||Il|,|12 



To introduce the R-charge of the BMN approach we single out two of the six scalars 
i = {1, . . . , 6}, which transform under the S0(6) R-symmetry group, say 05 and 
06; and form the complex combination 

Z(x) = -^(05(a:) + ^06(x)) (1.2) 

Then we define the R-charge, J, as the quantum number conjugate to the phase of Z. As 
mentioned above, operators which survive the BMN limit are characterized by having 
J very large and A — J finite. In practice this means that such operators contain a large 
number of Z-fields and a finite number of impurities in the form of fields not carrying 
R-charge such as 0i, 02, 03 and 04. In A/" = 4 SYM and in particular in its BMN sector 



■^It appears that this hmit is the same for gauge groups SU(N) and U(N). In this paper we shall 
be considering gauge group U(N). 



1 



the space-time dependence of two- and three-point functions is fixed by conformal 
invariance. At the classical level the calculation of such correlators then reduces to 
the calculation of expectation values in a zero dimensional Gaussian complex matrix 
model. For protected operators this statement trivially remains true when interactions 
are included and for non-protected operators a similar simplification can be obtained 
even at the quantum level if one introduces effective vertices . (This procedure has 
so far only been implemented at one- loop level.) In reference fl^ it was proposed that 



only two-point functions of appropriately defined multi-trace operators would have 
a string theory interpretation and this point of view has been supported by gauge 
theory calculations 0,0. This implies that extracting information about pp-wave 



strings from the gauge-theory reduces to determining the expectation value of traces 
of words in a zero-dimensional Gaussian complex matrix model. So far genuine matrix 
model techniques have only been exploited in the calculation of a very limited set of 
expectation values whereas the major part of those obtained were determined 
by combinatorial means. For higher genera combinatorial arguments become very 
involved. From the string theory point of view higher genera contributions are most 
interesting because they encode information about string interactions. So far gauge 
theory calculations were only pursued up to and including genus one. 

In the present paper we shall develop a loop-equation based, analytic strategy 
which allows us to calculate by recursion expectation values of products of arbitrary 
traces of words in a Gaussian complex one-matrix model to any order in the genus 
expansion. The outline of our paper is the following. First, in section 2 we explain in 
more detail how to reduce the calculation of two-point functions in A/" = 4 SYM to the 
calculation of matrix model expectation values, focusing on the two-point function of 
the so-called BMN operators. Next, in section 3 we introduce the notation necessary 
for our matrix model investigations and list the matrix model expectation values which 
are needed to find respectively the genus one and the genus two, one-loop correction to 
the anomalous dimension of the BMN-operators. In section 4 we derive the two basic 
relations on which all our considerations are based; the split and merge rule respectively. 
As a first application of these rules we reproduce in section 5 all the matrix model 
expectation values needed for the above mentioned genus one calculation by purely 
analytic computations. Subsequently, in section 6 we determine the correlators needed 
for the genus 2 calculation and finally in section 7 we show how our strategy allows 
us to find the expectation value of traces of arbitrary words to any order in the genus 
expansion. Section 8 is devoted to correlators which can be calculated exactly and 
section 9 contains our conclusions. 

Note: As we were completing our manuscript a related, interesting paper appeared 
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where another loop equation based technique is apphed to the study of the (planar) 
BMN-limit |6|. 



2 From J\f = ^ SYM to matrix model 

The field content of A/" = 4 SYM in four dimensions consists of the scalars 0i(a;), 
i G {1, . . . , 6}, a space-time vector A^{x) and a sixteen component spinor ip{x). These 
fields are Hermitian N x N matrices and can be expanded in terms of the generators 
of the gauge group U{N), for instance 

{<P^)ap{x) = J2 ^t{x)T:p (2.3) 



The generators are normalized as follows 



a=0 



tr [T\ T'] = ^^P^is = SsJp-r (2-4) 

a=0 

and the Euclidean action reads 



+^V^r^D^^-^V^r,[0„^]^ (2.5) 

where F^^^y = d^A^ — di^A^ — i[A^,Ay] and the covariant derivative is -D^0i = 
df_i4>i — i[Afj_,(j)i]. Furthermore, (F^jFj) are the ten-dimensional Dirac matrices in the 
Majorana-Weyl representation. Working in Feynman gauge, the propagators of the 
scalar fields take the form 

9ym 

{{(Pi)al3{x){(t)j)y5{0)) = -^^^^tj^asSfi^ (2.6) 

and in particular (cf. eqn. ( p..2|) ) 



{Z^pix)Z^sm = 5as5p^ (2.7) 

{Z^fsix)Z^sm = {Z^fsix)Z^sm = (2.8) 

Operators 0{x) which belong to the BMN sector of A/" = 4 SYM are characterized by 
containing a large number of Z-fields and a finite number of impurities in the form of 
fields not carrying R-charge. As an example, let us consider the most studied, so-called 
BMN-operator 

O'uA^) = ^^i=f^Ee^^^^"/'tr {Mx)Z'{x)Mx)Z'~'{x)) (2.9) 
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From the Feynman rules ( p.6[ ), ( ^.T] ) and ( ^.8| ), (or alternatively from conformal invari- 
ance) it follows that the tree level two-point function of BMN-operators can be written 
as 

{0(2,ni^)0i,^^m = [^) E e^-^(")/^(tr(0iZP02^^-P)tr(0iZ^-''02^^)) 

^ / p,q=0 

where the space-time independent matrix valued fields, (pi and Z should be contracted 
using the following Feynman rules 

((0j)a/3(0i)75) = SijSasSfj-y (2-10) 

(ZafsZ-ys) = SasS/s-r (2-11) 
The contraction of the 0-fields can easily be done by hand and we are left with 

{0(,^MO(,^mm = Yl e^-^("^-'"^)/^(tr (Z^-^Z^-^)tr(Z^Z^)) (2.12) 

V / p,q=0 

Now, the remaining expectation value can be identified as an expectation value in a 
zero dimensional Gaussian complex matrix model, namely 

{ti {Z -P Z ^'') ti {Z^Z'^)) = j dZdZ exp {-ti (ZZ)) ti {Z-^-^Z-^"') ti {Z^Z'^) (2.13) 
Here dZdZ is defined as 

dZdZ = n ^R^^^I"^^ (2.14) 
J. J. ^ 

such that / dZdZe-^'^^^^ = 1. The identification (|1|) holds because the matrix 
model measure ( |2.14{ ) combined with the Gaussian action precisely give rise to the 



contraction rule ( p.ll| ). The action and the measure carry a U{N) x U{N) symmetry 
corresponding to the transformation Z UZV'^ with U and V unitary. Expectation 
values of operators likewise carrying this symmetry, i.e. traces of products of {ZZ) 
can be calculated even for arbitrary U{N) x U{N) invariant potential order by order 
in the genus expansion using loop equations |T7|. Expectation values of operators 



consisting of products of traces involving only Z's or Z's can likewise be obtained by 



well-established methods, namely by character expansion |jl8|,[T^ or by the method of 



Ginibre ||2^. Notice that the object appearing in (p.l3|) does not belong to either of 



these classes of correlators. The aim of the present paper is to develop a method which 
allows us to deal with general correlators composed of traces of arbitrary words of Z 
and Z. 

It is obvious that any tree-level two point function of operators in the BMN sector of 
= 4 SYM can be reduced in the above manner. One pulls out the space-time factor, 
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contract by hand the finite number of impurities and one is left with a matrix model 
expectation value. By making use of so-called effective vertices one can also reduce 
one-loop corrections to two-point functions to matrix model expectation values 0,^. 

As explained in the introduction correlation functions in the BMN sector of A/" = 4 
SYM have an expansion in powers of j^, the genus counting parameter, and we are 
interested in determining (at least) the first terms in this expansion. For that purpose 
it is convenient to decompose our matrix model expectation values into connected and 
disconnected parts, f. inst. 

( tiiZPZ") triZ^-PZ^-")) = ( tr (Z^'Z^)) ( tr (Z-^-^'Z^-'?)) + ( tr (Z^'Z'?) ti {Z^-PZ-^-'')),onn 

as the connected part is down by a factor of compared to the disconnected one. 

Furthermore, it is convenient to work with generating functionals for expectation 
values in stead of working with the expectation values themselves. For instance, let us 
define 

W^i,i(xi,yi;x2,2/2) = (tr ( \ ^ \ ) tr ( \ _ )) (2.15) 

\ \xi-Zyi-ZJ \X2- Zy2- Z J / ^^^^ 

where xi, yi, X2, y2 are to be viewed as auxiliary variables. Then we have 

, _ — iTrn —iivm iivn iirm , 

Wi^i{Xe—,Xe^;Xe—,Xe—)= (2.16) 

oo J 

^in(m-n) ^(XX)"'^"^ ^ ( tr (Z^-^Z^^^) tr {Z'P Z")) ^^^^e^'""^"'^-""^^/-^ 
j=o p,q=o 



which immediately allows us to extract the sum appearing in (|2.12|) also known as 



the tree-level mixing matrix. So far, the tree-level mixing matrix has been calculated 
to order (genus one) in |]^,|^ and to order (^j^^ in 0- Furthermore, the one- 
loop correction to the two-point function was calculated to genus one in In the 
following section we list the matrix model expectation values or rather the generating 
functions needed for that computation. We hkewise list the ones needed to extend that 
calculation to genus two. Later we shall determine all of these functions.. 

3 Definitions and Notation 

We consider a complex Gaussian matrix model whose partition function is given by 

Z = [ dfxe-^ = [ dZdZ e-^*--^^ (3.17) 



where the integration runs over complex N x N matrices. Note that there appears a 
factor of N in front of the action. This factor is introduced only for convenience and 
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can easily be scaled away in the final results. Let us define the following generating 
functionals, also denoted as loop functions. 



1 

X 



-/ N 1 / 1 

Lij[y) = — ( tr =■ 

\ y-Z 



1 

y 



(3.18) 

(3.19) 
(3.20) 
(3.21) 
(3.22) 

y I conn 

We have normalized these functions so that their leading term in the large-N expansion 
is of order one. Knowing the leading order contributions to these functions for large N as 
well as the next to leading order contribution to Wi{x, y) suffices for the calculation of 
the one-loop, genus one correction to the anomalous dimension of the BMN operators. 
However, we shall be interested in more general loop functions. We define 



1 / 1 1 

Wi{x, y) = — (ti = 

A'\ X — Z y — Z 



/ A 1 / 1 1 1 1 

W2{x, y,x ,y) = — (ix ^- ^ 

N \ X — Zy — Zx— Zy' — Z 



■tr 



W^i,i(x,y;x,y ) = ( tr 
Ui{x;x',y') 



x-Zy-Z x' -Zy' -Z 

1 11 
tr — tr 



• • ■ ,Xl,l^,yi,h'^ ■ ■ ■■,Xn,l,yn,l, ■ ■ ■ ,Xn,l„,yn,l„) 



(3.23) 



\j=i 



2 ^ yj,i 




l^^/i?..,in(^l,l'^l,l' • • ■ ^^IM^yhh'^ ■ ■ ■ Xn,l, yn,l, ■ ■ ■ ,Xn,ln,yn,ln) 



]\f2h 
h=0 

This function is invariant under permutation of the various traces, under cyclic per- 
mutation of the factors inside a given trace and it is changed to its complex conjugate 
under x ^ y. We can represent it with a Young diagram like graph as follows with 

/l < /2 • • • < 



X 

1,1 






X 

1,1, 




X 

2,1 






X 

2,1, 


^2,1, 



X 

n,l 






X 

n,I 

n 


y 
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We also define 

Ui^,...,i„{x] Xi^i, yi^i, . . . , Xi^ij, . . . ; Xn,i, yn,i, ■ ■ ■ , 
= AT"-! 



(3.24) 



tr 



X — Z 



JpK^iu-Jn^^'' 1,1/1,1, . . . , xi,/j, ?/i,/^; . . . ; a;„,i,?/„,i, . . . ,Xn,i„, Vnu) 



h=0 



which we can similarly represent with a Young diagram like graph 



X 






X 

1,1 


^1,1 




X 

1,1, 






X 

2,1 


\l 




X 

2,1, 




• • 

• 

• 




X 

n,l 






X 

n,l 

n 


y 



where h < h ■ ■ ■ ^In- Calculating the one-loop, genus two correction to the two-point 
function of BMN-operators would require the knowledge of the third order contri- 
bution to Wi{x,y), the next to leading order contribution to the functions (|3.2CI|) - 
( ^.22| ) as well as the leading order contribution to the functions U2{x; xi,yi,X2,y2), 
W3{xi,yi,X2,y2,X3,y3), Ui^i{x;xi,yi;x[,y[) and Wi^2{x,y;xi,yi,X2,y2)- In section 6 
we shall show how to determine these and in section 7 we shall describe a general 
strategy for determining any multi-loop function or equivalently any expectation value 
of traces of words to any order in the genus expansion. 



4 Loop Equations 

All our computations will be based on two simple rules which can be derived by loop 
equation techniques, based on the fact that the matrix model partition function is 

Here we restrict ourselves to considering the 



invariant under field redefinitions 22 



case of a complex matrix model with a Gaussian potential which is the case of interest 
for the BMN sector of A/" = 4 SYM. This is, however, just an almost trivial application 
of a method which works under much more general circumstances and will be presented 
for a Hermitian two- matrix model with arbitrary U{N) invariant potentials in [p^] . 
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4.1 Split Rule 

Consider the field redefinition 



Z ^ Z + eA^—B (4.25) 

X — Z 



This redefinition gives rise to the following change of the measure 

S (dZdZ) = 2Re ( e tr A^— tr -^—B] dZdZ + 0(e^) (4.26) 
\ X — Z X — Z J 

If A or 5 depends on Z ot Z there will be additional contributions which are obtained by 
applying the usual chain rule in combination with the split and merge rules. Obviously, 
under ( [4.25| ) the action changes as 



6S = 2NRe [ e tr (^ZA-^B ] ] (4.27) 



The relations (|4.26|) and ( |4.27|) hold for arbitrary complex e, in particular for e purely 



real or purely imaginary. Therefore we conclude 



tTA^—tT^—B)=N(tTZA^—B) (4.28) 
X — Z X — Z \ X — Z 



4.2 Merge Rule 

Here we consider the following field redefinition 



Z ^ Z + tr — — B (4.29) 

X — Z 



for which the change in the measure is 



(5 [dZdZ) = 2Re [etrA^—B^—] dZdZ + 0{e'^) (4.30) 

y X — Z X — Z J 

Again, if A or i? depends on Z or Z there will be additional contributions which are 
obtained by applying the usual chain rule in combination with the split and merge 
rules. The change of the action is obvious and our final merge rule reads 

tr ( A^—B^—] ) = N ( tr ZA tr -^—B) (4.31) 
\ x — Z x — Z) I \ x — Z I 

5 Functionals needed for the one-loop, genus one 
computation 

In this section we shall determine the leading order contribution for large N to the 
loop-functions (|3.19| )- (|3.22| ) as well as the next to leading order contribution to ( |3.19|) . 
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As mentioned above these are the objects needed for the computation of the one-loop, 
genus one correction to the anomalous dimension of the BMN-operators. In section 8 
we will show that the functionals ( |3.19 ) and (|3.22 ) can be calculated exactly but to 
expose the completeness of our loop equation method we shall derive the leading order 
contributions to these below as well. 



5.1 u{x) and Wi{x,y) to leading order 

Considering the field redefinition 6Z = we easily get 



= — 1 + xuj(x) 



which is true to all orders in and gives 



X 



(5.32) 



(5.33) 



This result of course trivially follows from symmetry arguments. Next, we make use of 



the field redefinition SZ 



1 1 

x—Z y—Z 



and obtain 



X y 



= y 



uj{x)Wi{x, y) + j^Ui{x] x, y) = yWi{x, y) - u{x) 



(5.34) 



Above, a space between two Young diagrams signifies multiplication of the correspond- 
ing functions. To leading order in we can neglect the second term in ( ^.34| ) and we 

get 

Wf'\x,y) 



xy — 1 

which is in accordance with the simple combinatorial result 

^{ti Z^Z^) = 1 + 0{1/N^) 

5.2 W2, Ui and Wi^i to leading order 

Performing the change of variable 6Z = — — 

o o xi—Z yi~Z X2 — Z y2—Z 



leads to: 



(5.35) 



(5.36) 



9 



1^2(3^1, 2/1, X2, 1/2) (2/2 - — ) 

Xi 



(5.37) 



:i + W^i(x2,2/2)) 



Wi{xi,yi) -Wi{x2,yi) 



X2 — Xi 

1 W^i,! (2^1,2/1; 3^2, 2/2) - Wi^i{x2,yi;x2,y2) 



+ 7?2^2(a;i;a;i,2/i,X2,2/2) + ,,2 

jS'^ X2 — Xi 

In the equation above we have used fractional decomposition to express the quantities 
represented by the two last Young diagrams in each line in terms of usual W functions 
with an even number of arguments. From equation (|5.37|) we can easily find the leading 
contribution to W^2(a^i, 2/i, a^2, 2/2) for large N, namely 



Wi°\xi,yi,x2,y2) 



xiyiX2y2 



(xi2/i - l)(x22/i - l)(a;i2/2 - l)(a;22/2 - 1) 
which reproduces the combinatorial result 
1 



tr Z-^-^'Z-^-'^Z^Z") = 1 + Min [p, q,J-p,J-q]+ 0(1/ 



(5.38) 



(5.39) 



Next, carrying out the change of variables SZ — xi-Zyi-z 
simple relation 



^ tr „ -"^ ^ ._ ^ ^ we get the following 



X 














X 

1 


y 

1 


1 

~ X 


X 

1 


y 

1 


X 


y 

1 



i.e. to leading order 



f/i(x;xi,2/i) = -W2ixi,yi,x,yi) 

X 



U[°\x;xi,yi) 



xiyl 



(xi2/i - lY{xyi - 1)2 
which is in agreement with the combinatorial result 

( tr Z^ tr Z'Z'-^)^^^^ = p( J _ p + 1) + 0(l/iV2) 



Finally, we consider the field redefinition bZ 
the following equation for W\;y{x\^ 2/ij 2:2, 2/2) 



1 1 tr 1 1 



(5.40) 
(5.41) 

(5.42) 
which leads to 
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X 

1 


y 

1 


X 

2 


y 

2 



+ 



X 


y 


1 


1 



X 

1 




X 

2 


y 

2 



+ 



X 

1 




X 

1 


y 

1 


X 

2 


y 

2 



+ 



X 


X 


y 


X 


y 


2 


1 


1 


2 


2 





X 

1 


y 

1 




X 

1 




= yi 


X 

2 


y 

2 




X 

2 


y 

2 



W^i,i(a;i,2/i;a:2,2/2)(|/i - — ) 

= (1 + Wi{xi,yi)) Ui{xi;x2,y2) + 



(5.43) 



W'2(a;i,i/i,X2,?/2) - W2ix2,yi,X2,y2) 



X2 — X\ 



Now, making use of (|5.35| ), ( ^.38| ) and (|5.41|) we get 



W^U (a;i,yi;x2,i/2) 



where we note that 



x\X2y\y2{^ + xiX2yiy2{xiy2 + X2yi - 3)) 
{xiyi - iy{xiy2 - lY{x2yi - lY{x2y2 - 1)^ 



(5.44) 



1 xiy2 Xiy2 

[l + xiX2yiy2{xiy2 + X2yi - 3)) = det \ X2yi 1 X2yi 

X2yi xiy2 1 



(5.45) 



From (|5.44|) we can easily get the genus one correction to the tree-level mixing matrix 
of BMN operators coming from connected diagrams, namely (cf. eqn (|2.16| )) 



J 

= ^ ( tr {Z^-PZ^-'^) tr {Z^Z'^))^^^^ ^2niinp-mg) / J 
p,q=0 

dr 



(5.46) 



2m 



-inn/ J 



-inm/ J . 



iTTm/J\ i-K(n—m) 



dr 



.J +3 1 



1 + [2r cos 



{n-\-m)iT 



-3])e 



i-K(n—m) 



2,TCi (r Q—iiT{n—m,)/J'^2^ri^ g— i7r(n+m)/ J^2 gi7r(ra— m)/ 2 gi7r(n+m)/J^2 

It is obvious that the analyticity structure of the integrand depends on the values of n 
and m, more precisely we have 

• n = m = 0: A pole of order 8 at r = 1 

• n = and m 7^ (or m = and n ^ 0): Two poles of order 4 at r = e^*'^'^/'^ (or 
at r = e^^^'^™/-^) 

• n = m^OoTn = —m 7^ 0: One pole of order 4 at r = 1 and 2 poles of order 2 
at r = e^''^*'^'^/'^ Notice that the residues are not the same in the two cases. 
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• \n\ 7^ \m\ and n ^ ^ m: Four poles of order two at r = e^""^*-" '"•'/"^ and 

J. _ g±i7r(n+m)/J 

Strictly speaking, the conditions on n and m are to be understood modulo J but we 
always consider m <^ J . The fact that the evaluation of the mixing matrix has to be 
split into 5 separate cases follows immediately in the generating functional picture and 
evaluating the contour integral (|5.46| ) we easily reproduce the result of reference i.e. 



M^) ^ MO) ^ t5 



up to terms of order J'^. 



j_ 

40 

3-27r2 
24n47r4 
21-27r^ 

48 n^TT-* 
9 



32n47r4 

2n^-3nm+2m^ 
8 'n? m? (n—m)^ tt'* 



77, = m = 

n 0,171 = 

n = m and n 7^ 

= —m and n 7^ 
In I 7^ |m| and n 7^ 7^ m 



(5.47) 



5.3 1^1(0:,?/) to next to leading order 



Inserting the genus expansion ( 3.23|) into ( ^.34 ) we can easily determine the genus one 
contribution to Wi. From ( ^.341) and (|5.41|) we get 



X 



xy 



U[°\x; x,y) 



x'^y'^ 



{xy - If 



(5.48) 



6 Functionals needed for the one-loop, genus two 
computation 



In this section we shall determine the third order contribution to (|3.19|) , the next to 
leading order contributions to the functionals ( |3.2(JD -( ^.21 ) as well as the leading order 
contribution to the functions Ui^i, Wi^2, U2 and W3. We shall start by the latter ones 
and work our way toward the first ones. 



6.1 Ws, U21 W2,i and f/i 1 to leading order 

Considering the field redefinition 5Z 



1111 



x\-Z yi-Z X2~Z y2-Z X3-Z ys-Z 



we obtain 
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W^{xi,yi,X2,y2,xz,y^){y-i - —) (6.49) 

Xi 

, ..J I .y^2{x\,y\,X'2,y2)-W2{xz,yx,x2,y-^ 
= (1 + vri(x3,y3)) 

X3 — X\ 

+ W^2(a;2,2/2,a:3,?/3) 

X2 - Xx 

1 W^i,2(a;i, xa, ?/2, 2:3, 1/3) - W^i,2(a;2, y\\ ^2, y2, x^, y^) 



N"^ X2 — Xi 

1 ^i,2(a;3, 1/3; xi,yi, X2, 1/2) - ^^1,2(3^3, y^; a^3, l/i, a^2, 1/2) 

X3 — xi 



+ -^U^ixi, xi,yi, X2, y2, X3, y^) 

From here we can immediately get the genus zero contribution to from our knowl- 
edge of the genus zero contribution to Wi and W2. The result reads 



wl^\xi,yi,X2,y2,X3,y:i) 



^^hjy^'y^ ^) \ x^wi 1 1 



1 xiy2 1 
1 

xm 



(6.50) 



Next, performing the change of variables 5Z = -r^tr ^ ry ^r^ leads to the 



following simple relation 



x~Z xi~Z y\—Z X2—Z y2—Z 



X 




X 

1 


y 

1 


X 

2 


y 

2 



]_ 

X 



X 


y 


X 


y 


X 


y 


1 


X 


y 


X 


y 


X 


y 


1 


1 




1 


2 


2 


+ X 


1 


1 


2 


2 




2 



t/2(x; xi, 2/1, X2, 2/2) = - (W'3(xi, 1/1, X, 2/1, X2, 2/2) + W^3(xi, 2/1, X2, 2/2, x, 2/2)) (6.51) 



and from which we easily get U^^ by inserting ( 5.50 ). Furthermore, choosing the field 



redefinition 5Z 



1 1 tr^ 1 1 



x\~Z yi — Z X2—Zy2—Zx-i—Zyj, — Z 



we find 



13 



X 

1 



X 

1 


y 

1 












X 

2 


y 

2 


X 

3 


y 

3 


+ 


X 

1 


y 

1 



X 

1 




X 

2 


y 

2 


X 

3 


y 

3 



+ - 



X y 
3 3 



+ 



X 


X 


y 


X 


y 


X 


y 


+ 


X 


y 


X 


X 


y 


X 


y 


2 


1 


1 


2 


2 


3 


3 


2 


2 


3 


1 


1 


3 


3 





X 

1 


y 

1 








X 

1 




= yi 


X 

2 


y 

2 


X 

3 


y 

3 




X 

2 


y 

2 


X 

3 


y 

3 



W^i,2(a;i,?/i;a;2,?/2,a;3,l/3)(yi - — ) 

Xi 

= (1 + Wi{xi, yi))U2ixi; X2, 2/2, xs, ys) 

^ W3{X2, yi, X2, 1/2, 3^3, Vs) - Wsl^i, l/i, X2, t/2, 3^3, Vs) 

Xi — X2 

W3{x3, 1/3, X2, 2/2, X3, yi) - 1^3(3:3, 1/3, X2, 2/2, Xi.yi) 



(6.52) 



+ 



Xl - X3 



1 



+ J^Ui^2{.Xi] Xl, 2/1; X2, 2/2, 3^3, 2/3) 



where Wi^2 is now expressed in terms of already known quantities. Furthermore, making 
use of the chan^ 
simple relation 



use of the change of variable 5Z = tr — tr — we obtain the following 



X 




X 

1 


y 

1 


X 

2 


y 

2 



]_ 

X 



X 

2 


y 

2 




X 

1 


y 

1 


X 


y 

1 



1 



+ 



X 

1 


y 

1 




X 

2 


y 

2 


X 


y 

2 



^\,\{x] Xl, 2/1; X2, 2/2) = ^ (M^i,2(a;2, 2/2; a;i, 2/1, x, 2/1) + W^i,2(a;i, 2/1; X2, 2/2, a;, 2/2)) (6.53) 
6.2 iy2, L^i and VFi 1 to next to leading order 



Inserting the genus expansion ( p.23| ) and (|3.24|) into the relevant loop equations ( |5.37|) , 
( |5.40D and ( ^.431) we can easily determine the genus one contribution to f/i, W2 and 
W\;y. From ( |5.37|) we find (using Mathematica) 



W^2^^(a;i,2/i,a;2,2/2) 

J , ttKO) 



(6.54) 



iyi('^(xi,2/i)-l^r(x2,2/i) 



x\y2 - 1 



(i+iyr(^2,2/2)) 



X2 - Xl 
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X2 - Xi 



Wi°l{xi, yi; X2, 1/2) - Wi:°l{x2, yi, X2, 2/2) , (o) 



X2 - Xi 



+ '{xi, Xi,yi,X2,y2) 



rPol28(a;i,?/i,X2,l/2) 



where Pold(xi, yi, X2, ?/2) denotes a polynomial of degree d in xi, X2, y2- If we carry 
out the discrete Fourier transform and take the large-J limit we get 



(6.55) 



p,q=0 



dr 
27ri 



f j_ 

240 



48 n^TT-* 

315-1207r^ n^^m-K^n^ 
7680 nfivre 

3(n'^+m.'^)+47r^(n'^m^+m^n^)-15(m'*n^+n'*m^)-8m.'*n-^7r^ 
48 n2 m? {m-nY {m+n)47r6 



n = m = 

n 7^ 0, m = 

|?T,| = |m| and n 7^ 

|n| 7^ |m| and n 7^ 7^ m 



up to terms of order J^. 
Furthermore, from (|5.40|) 



Ijf'ix-.xx.yx) = -W^''{x,yi,xi,yi) 



(1), 



X 



(6.56) 



xiyl 



{1 + 4?/i(xi + x) — 34?/^a;xi 



(l-x^i)6(l-xiyi)6 
18|/^(a;xi(a; + Xi)) + y^xxi{—10x^ + 21xxi — lOx^) 

+2ylxxi{x^ — 6x^x1 — 6xxl + xf) + |/^x^x^(3x^ + 2xxi + 3x^)} 



Finally, from ( 5.45 ) one gets 
Wi:}{xi,yi;x2,y2) = 



(6.57) 



Xl 



(1 + Wl'\xx,yx))ui'\xi; X2, 1/2) + H^f ^(xi, yi)ui'\x,; X2, 1/2) 



- 1 

Wi^\xi,yi,x2,y2) -W^''{x2,yi,X2,y2) , rr(o) 



X2 — Xl 



+ Ul{{xi;xi,yi;x2,y2) 



I{^Axiy, - ly 



rPol36(xi,?/i,X2,l/2) 



In this case carrying out the discrete Fourier transform and taking the large-J limit 
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gives 



73 



181440 

7 _ 
lOn^TrS 

143 
2048n87r8 
7 



13 



107 



+ 



1 



768n67r6 
245 



320n^TT* 



480n27r2^ 



+ 



160n47r4 
13 



2016ri27r2 



4096n87r8 _^30 727i67r6 ' 1280ri47r4 
, 192m4n4(m-n)6(m+n)47r8 ^ Poli4(m, 77,) 



n = m = 

n 7^ and m = 

n = m and n 7^ 

n = — m and n 7^ 

|n| 7^ \m\ and n 7^ 7^ m 



where 

Poli4(m, n) 



+ Q3m^ + 9m^ n'^ - 13m^ - mn^ + 2n^^) 

- (m^-ra^)^ (12 - 18 n - 20 m^ra2 + 37 m^ri^ 

+38 + 37 it' - 20 ra^ - 18 m + 12 n^) tt^ 

+ m^n^ (m^ — n^)^ (2m^ — 3m?2 + 2?2^) vr^ (6.58) 

and where we have neglected terms of order J^. Notice that the expressions above 
constitute the contribution to the mixing matrix coming from connected diagrams only. 
If one includes also disconnected ones one reproduces the expressions given in 0.0 

6.3 Wi{x^y) to third order 

Making use of ( ^.34 ) and ( 6.56 ) we get 



X 



xy — 1 



U[^\x;x,y) 



2„-,2^ 



x^y^^l + I2xy + Sx^y 
{xy - 1)9 



(6.59) 



7 The general case 



From the examples above it should be clear how to choose the appropriate field redefini- 
tions needed for the derivation of the loop equation associated with a given generating 
functional. Here we shall write down the most general loop equations and show how 
they allow us to determine recursively any multi-loop correlator, i.e. any expectation 
value of traces of words to any order in the genus expansion. 



7.1 One-trace functions 

Considering the change of variable 5Z 
the following relation 



1 



1 



xi-Z yi-Z ' ■ ■ xi-Z yi-Z 



~ with / > 2 we obtain 



[l/xi - yi)Wi{xi, 



yi) 



(7.60) 



detailed comparison is not possible in the general |7i| 7^ \m\ case and there seems to be a 
discrepancy in the two last terms of the n = —m case. 
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^ Wk-i{xi,yi, . . .,yk-i) - Wk-i{xk,yi, . . .,yk-i] 



Xi — Xk 
k=2 ^ 

i5k,i + Wi-k+iixk, yk,---, xu yi)) + 0{j^) 

Here Wi is expressed entirely in terms of Wk with k < I — 1. Having determined the 
planar contribution to Wi{x^y) (cf. eqn. ( [5.35D ) we can by means of ( |7.6(J| ) determine 



recursively the planar contribution to any one-trace function Wi. From the structure 
of the recursion relation and the explicit expression for wf\x, y) it follows that planar 
one-trace functions only have singularities in the form of single poles. More precisely 
we have 

yi, . . . , xz, yi) = — JIi^i^Pol(a:i, yi, . . . , xu m) (7.61) 

where Pol(xi, . . . ,xi, yi) is a polynomial of degree / — 2 in each of its variables. In 
the case / = 3 (and trivially in the case / = 2) this polynomial could be expressed 
as a determinant (cf. equations ( ^.38|) and ( |6.50| )) but it does not seem that a similar 



simplification occurs for higher values of /. It would be most interesting, though, to 
find a closed expression for Wi for general /. 

7.2 Multi-trace functions 

In the case of the VT-functions it is convenient to consider separately the cases h = 1 
and Zi > 1. For h = 1 and (n > 2) we have 

(yi,i —)Wi^i^^„j„{xi^i,yix,X2,i,...,yn,iJ = (7.62) 

Xl,l 

(1 + W{Xi^i, yi,l))t/i2,...,/„(Xi,i; X2,l, . . . , yn,ij 
n Ik ^ 

, yk,j-i, 1/1,1, Xk,j, ykj, ■ ■ ■) 

k=2 j=l ^^^^ ^I'l 
— W^h,-,'fc+l, ■■■,'.!(■ ■ ■ 5 VkJ-l, Xk,j, 1/1,1, Xkj, yk,j, ■ ■ •)} 

+-^t^i,«2,...,«n(a;i,i;a;i,i,i/i,i;x2,i, . . .,?/„,«„) 
whereas for li> 1 the relevant loop equation reads 

{yi,h — ^Wi^,..u{^i,i,-,ynu) (7-63) 

a;i,i 



h 

^ iy,i+i_j,..,/„(a;i,j, .., ?/„,iJ(l - 5„,ij 
Wj_i(xi,i, 1/1,1, .., i/i,j-i) - iyj_i(xi,j, 1/1,1, .., yij-i) 



h 



Xl^j — Xii 



i=2 
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k+h,..,lnix2A, ..,|/fcJ-l,Xi,i, ..,yi^l^,Xk,j, ■■,yn,u) 

k=2 j=l ^^'J' ^1-1 
-Wll,..,lk+ll,-,lni^2,l, ■■, a^fcj, .., yi,h^^k,j, ■■,yn,in)} 

+lv2 Z {^j-i,Zi+i-i,i2,..,Zn(a;i,i, 1/1,1, ..,yi,j;xij, ..,?/„,/„) 

■ o Xl j Xi i 

-W^i-iA+i-i,..,«n(a;i,i, yi,i, Z/ij-; .., i/n,i„)} 
+^t/zi.../„(a;i,i;a;i,i,... 

These equations are to be supplemented by the loop equations for the ^/-functions 
which take the simpler form 

X Ui^,.. j„{x; xi^i,. .. ,yn,ij (7.64) 

= ^ ^ W^/i,...,ife+i,...,/„(a;i,i, • • • ; Xk,i, • • . , yk,j, x, ykj, Xkj+i, • • • , 2/^,/^; • • • , yn,i„) 
k=i j=i 



The relations ( [7.62| ), (|7.63 ) and (|7.64|) constitute a triangular set of equations which 



allows us to determine any multi-loop function to any order in the genus expansion. 
For a finite number of loops and finite genus only a finite number of operations are 
needed. Below we shall make this statement more precise. 

First, let us introduce an ordering of the Young diagrams representing ly-functions. 
A Young diagrams is said to be smaller than a Young diagram Yi^^ (repre- 

senting loop-functions and respectively) if: 

m n 

m + ^ki <n + ^li (7.65) 

i=l 1=1 



Next, let us consider the loop equations ( [7.62|) and (|7.63|) (using ( [7.64|) ): it is clear that 



all the leading order diagrams on the RHS are smaller than the diagram on the LHS. 
This means that at the planar level a ly-function corresponding to a certain Young 
diagram can be expressed entirely in terms of planar l^-functions corresponding to 
smaller Young diagrams. We thus get a closed equation for any genus zero ly-function 
and clearly also for any genus zero [/-function (cf. equation ( [7.64] )). 

Proceeding to higher genera, we have in our loop equations (|7.62|) and ( [7.63| ) two 
types of terms which carry a factor l/N"^; [/-terms and IV-terms. Compared to the 
object we are interested in, the l^-terms correspond to Young diagrams where one 
extra line has been added while the number of boxes has been kept fixed. The U- 
terms, on the other hand, correspond via ( |7.64| ) to Young diagrams where two extra 



boxes have been added while the number of lines has been kept fixed. This means that 
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the genus g contribution to a W^-function described by a Young diagram with 2K boxes 
and n hues can be expressed entirely in terms of genus zero W^-functions corresponding 
to Young diagrams having at most 2K + 2g boxes and n + g — 1 hues. Clearly, we 
have a triangular set of equations which allows us to determine any expectation value 
of traces of words to any order in the genus expansion. 



8 Exactly calculable correlators 



As shown in reference [0,0] it is possible to find exact expressions for the expectation 



values encoded in the following generating functionals 

= AT^-i/tr^-... tr^-tr-1^ ) (8.66) 



as well as [14 



G(a;i, xs, yi, ya) = A^^ ( tr tr tr tr ) (8.67) 

xi - Z X2- Z yi- Z y2- 



conn 



From the generating functionals ( |8.66| ) and ( |8.67| ) it is possible using again loop equa- 



tions to derive exact expressions for yet other generating functionals. The functions 
Wi{x, y) and Ui{x'\ x, y) can be determined in full generality whereas the method only 
gives the remaining W- and [/-functions in certain limits where typically a number 
of their arguments are sent to oo. To obtain the all genus version of Wi{x,y) one 
considers the field redefinition 6Z = tr which leads to the following first order 

x—Z y—Z ° 

differential equation 

-dyWx{x,y)=xEx{x,y) (8.68) 
which is to be supplemented by the boundary condition 

Vri(x, ?/) — > — , as Ixl, lyl ^ oo (8.69) 
xy 

The equation ( ^.68 ) is of course nothing but the generating functional version of the 



simple relation 

{iiZHiZ^) = J {iiZ^-^Z^-^) (8.70) 

which implies 

wi 1 1 r(Ar + j + i)! m \ 

'^^'^^ A^^x-^+V+M-^ + l)(-^ + 2) I (AT-l)! (AT- J- 2)! J 



j=o 

which can also be written 



m(x,y) = 5^-^/2fc+i(xy) (8.72) 



/c=0 
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with /i(a;) = l/(a; — 1) and: 



„ 1 d d^ 



dx^ '^'^^^ 1 — X dx dx'^ ^ ^ 

and fk{x) = 0{l/x^) for large x. The coefficient of l/(a; — 1)^*"'^^ as x ^ 1 is {2k — 3)\\/k. 

In the case of Ui{x'; x, y) one chooses the field redefinition 5Z = tr ^j^^ tr 
and obtains 

- dyUi{x';x,y) = xH2{x',x,y) (8.74) 
and the appropriate boundary condition in this case reads 

Ui{x';x,y) \ as |?/| ^ cx) (8.75) 
[X yj'^x 

Expressed in terms of expectation values ( |8.74| ) reads 



( tr tr Z^ tr Z^+^) = {J + K) {tr Z^ tr Z^+^-^Z^-^) (8.76) 
which has the obvious generalization with J = X]i=i 



tr 



Z-^YltiZ-^A = J I tr Z-^-^Z-^'-^ Yl tr Z-^^ \ (8.77) 



1=1 



9 Conclusion 

With this work we have added pp-wave physics and A/" = 4 SYM to the long list of 
areas where classical matrix model techniques have proven very efficient. 

As explained in section 2, evaluating a typical correlation function in the BMN 
sector of A/" = 4 SYM can be reduced to evaluating the expectation value of a product of 
traces of words in a zero- dimensional Gaussian complex matrix model and subsequently 
carrying out a discrete Fourier transformation. Simple correlators can be obtained by 
purely combinatorial arguments but such arguments become more and more involved 
(and correspondingly less and less reliable) the more words enter the correlators and 
the higher genus one is aiming at. With our loop equation based technique, however, 
we can by analytical manipulations reach any multi-word correlator to any order in the 
genus expansion. Furthermore, by working with generating functionals we trade the 
process of Fourier transformation for simple contour integration. 

There are several directions of investigation where our technique would be most 
useful. One is the investigation of operators with more impurities than the traditionally 
studied BMN operators of equation ( p.9|) . Such operators would correspond to string 
states with many oscillators excited and determining their correlators would imply 
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evaluating expectation values of many letter words. Such words are encoded in U- and 
W^-functions whose sub-scripts are large and these are of course accessible with our 
method. 

As mentioned several times our method would also allow us to calculate the higher 
genera, one- loop corrections to the anomalous dimension of the BMN operators ( |2.9| ). 
In fact, we have already evaluated all expectation values needed for the genus two 
calculation. As pointed out in [l^, completing this calculation would allow one to check 
whether the effective string coupling constant in the pp-wave/BMN correspondence is 
indeed 5^2 as suggested in Finally, it is possible that like in the one- loop case 

two- and higher loop computations on the gauge theory side can be reduced to pure 
matrix model computations and then obviously our method will again be in demand. 
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